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ABSTRACT 

Magic labelings of graphs are studied in great detail by Stanley in [18| and [19|, and Stewart in [20| and 
ITfl . In this article, we construct and enumerate magic labelings of graphs using Hilbert bases of polyhedral 
cones and Ehrhart quasi-polynomials of polytopes. We define polytopes of magic labelings of graphs and 
digraphs. We give a description of the faces of the Birkhoff polytope as polytopes of magic labelings of 
digraphs. 
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Figure 1: A magic labeling of Kq [21 1. 



1 Introduction. 

Let G be a finite graph. A labeling of G is an as- 
signment of a nonnegative integer to each edge of 
G. A magic labeling of magic sum r of G is a la- 
beling such that for each vertex v of G the sum of 
the labels of all edges incident to v is the magic 
sum r (loops are counted as incident only once) 

OUGllEOlED- Gra P hs with a ma § ic labeling are 
also called magic graphs (see Figure ^for an ex- 
ample of a magic labeling of the complete graph 
Kq of magic sum 40). 



We define a magic labeling of a digraph D 
of magic sum r to be an assignment of a nonneg- 
ative integer to each edge of D, such that for each 
vertex v% of D, the sum of the labels of all edges 
with Vi as the initial vertex is r, and the sum of 
the labels of all edges with Vi as the terminal ver- 
tex is also r. Thus magic labelings of a digraph is 
a network flow, where the flow into and out of ev- 
ery vertex, is the magic sum of the labeling. If we 
consider the labels of the edges of G as variables, 
the defining magic sum conditions are simply lin- 
ear equations, and the set of magic labelings of G 
becomes the set of integral points inside a pointed 
polyhedral cone Cq lfT31 . Henceforth, we call Cq 
the cone of magic labelings of G. A Hilbert basis 
of Cq has the property that any magic labeling of 
G can be expressed as a linear combination with 
nonnegative integer coefficients of the elements of 
the Hilbert basis (see Q, 0, or fl5l ). An irre- 
ducible magic labeling of a graph is a magic la- 
beling that cannot be written as a sum of two other 
magic labelings. The minimal Hilbert basis of Cq 
is the set of all irreducible magic labelings and 
is unique [15]. Henceforth, when we refer to the 
Hilbert basis in this article, we mean the minimal 
Hilbert basis. 

Let v\,V2, ■ ■ ■ ,v n denote the vertices of G 
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and let e^, ej 2 , . . . , ei m . denote the edges of G 
that are incident to the vertex vi of G. Consider 
the polytope 

rn, 

V G = {LeC G C R", ^L(e 4j ) = l-i = 1, . . . ,n}. 

3=1 

We will refer to Vq as the polytope of magic 
labelings of G and denote Ha(r) to be the num- 
ber of magic labelings of G of magic sum r. Then, 
Hc{r) is the Ehrhart quasi-polynomial of Vg IT6l . 

We define a polytope of magic labelings 
of a digraph D as follows. Let , . . . , ej m . de- 
note the edges of D that have the vertex vi as the 
initial vertex and let , /j 2 , . . . , /j s . denote the 
edges of D for which the vertex Vi is the terminal 
vertex, then: 

v D = {l g c D c R9, E, m i^K) 
= E^i^(4) = 1;^ = i,---, n}. 

Let Hi)(r) denote the number of magic la- 
belings of D of magic sum r. Then like before, 
-ffo (r) is the Ehrhart quasi-polynomial of Vd- We 
now connect the magic labelings of digraphs to 
magic labelings of bipartite graphs. 

Lemma 1.1. For every digraph D, there is a bi- 
partite graph Gd such that the magic labelings 
of D are in one-to-one correspondence with the 
magic labelings of Go- Moreover, the magic sums 
of the corresponding magic labelings of D and 
Gd are also the same. 

Proof. Denote a directed edge of a digraph D 
with Vi as the initial vertex, and Vj as the ter- 
minal vertex, by eij. Let L be a magic label- 
ing of D of magic sum r. Consider a bipartite 
graph Gd in 2n vertices, where the vertices are 
partitioned into two sets A = {ai, . . . , a n } and 
B = . . . , b n }, such that there is an edge be- 
tween a-i and bj, if and only if, there is an edge 
in D. Consider a labeling Lq d of Gd such that 
the edge between the vertices ai and bj is labeled 
with L(eij). Observe that the sum of the labels 
of the edges incident to a-i is the same as the sum 
of the labels of incoming edges at the vertex Vi of 




Figure 2: A digraph Dp and its corresponding bi- 
partite graph Gd q - 



D. Also, the sum of the labels of edges at a ver- 
tex bj is the sum of the labels of outgoing edges 
at the vertex Vj of D. Since L is a magic label- 
ing, it follows that Lq d is a magic labeling of 
with magic sum r. Going back-wards, consider a 
magic labeling L' of Gd- We label every edge 
of D with the label of the edge between and bj 
of Gd to get a magic labeling Ld of D. Observe 
that L' and L d have the same magic sum. Hence, 
there is a one-to-one correspondence between the 
magic labelings of D and the magic labelings of 
Gd- D 

For example, the magic labelings of the Oc- 
tahedral digraph with the given orientation Do in 
Figure |2] are in one-to-one correspondence with 
the magic labelings of the bipartite graph Gd q - 

A graph G is called a positive graph if for 
any edge e of G there is a magic labeling L of G 
for which L(e) > ED. Since edges of G that 
are always labeled zero for any magic labeling of 
G may be ignored to study magic labelings, we 
will concentrate on positive graphs in general. We 
use the following results by Stanley from [ 18 1 and 
fT9l to prove Theorems 1 1 . 41 and 1 1 . 5 1 and Corollary 

Em 

Theorem 1.1 (Theorem 1.1, |19|). Let G be a 

finite positive graph. Then either Hc(r) is the 
Kronecker delta 5q t or else there exist polynomi- 
als Ic{ r ) an d JG( r ) such that Hc(r) = Icif) + 
{-l) r J G {r)forallr € N. 

Theorem 1.2 (Theorem 1.2, |19|). Let G be a 

finite positive graph with at least one edge. The 
degree ofHc{r) is q—n+b, where q is the number 
of edges of G, n is the number of vertices, and b is 



the number of connected components of G which 
are bipartite. 

Theorem 1.3 (Theorem 1.2, |18|). Let G be a 

finite positive bipartite graph with at least one 
edge, then Hg(t) is a polynomial. 

We now conclude that Hf)(r) is a polyno- 
mial for every digraph D. 

Corollary 1.3.1. If D is a digraph, then Hr>(r) is 
a polynomial of degree q — 2n + b, where q is the 
number of edges ofD, n is the number of vertices, 
and b is the number of connected components of 
the bipartite graph Gd- 




Figure 3: The edge graph of Vr 3 - 



Proof. The one-to one correspondence between 
the magic labelings of D and the magic labelings 
of Gd, implies by Theorem 1 1.31 that Hd{t) is a 
polynomial, and by Theorem 11.21 that the degree 
of Hf)(r) is q — 2n + b, where b is the number of 
connected components of Gd that are bipartite. □ 

Consider the polytope V := {x\Ax < b}. 
Let c be a nonzero vector, and let 5 = max Ac < 
b}. The affine hyperplane {x\cx = 5} is called a 
supporting hyperplane of V. A subset F of V is 
called a face of V if F = V or if F is the inter- 
section of V with a supporting hyperplane of V. 
Alternatively, F is a face of V if and only if F is 
nonempty and 

F = {x e V\A'x = b'} 

for some subsystem A'x < b' of Ax < b. See lU31 
for basic definitions with regards to polytopes. 

Therefore, a face of Vg is a polytope of the 

form 




7! ffl m m if 




vi m m m 




Figure 4: The face poset of Vt 3 ■ 



ofG, n is the number of vertices, and b is the num- 
ber of connected components of G that are bipar- 
tite. The d-dimensional faces of Vg are the d- 
dimensional polytopes of magic labelings of pos- 
itive subgraphs of G with n vertices and at most 
n — b + d edges. 



{L eV G ,L{e ik ) =0;e ik e E }, 

where E$ = {e^ , . . . , ei r } is a subset of the 
set of edges of G. 

Theorem 1.4. Let G be a finite positive graph 
with at least one edge. Then the polytope of magic 
labelings of G, Vg is a rational polytope with di- 
mension q — n + b, where q is the number of edges 



Observe from Theorem 11.41 that there is an 
edge between two vertices V{ and Vj of Vg if and 
only if there is a graph with at most n — b + 1 
edges, with magic labelings V{ and vj. The edge 
graph of Vr 3 is given in Figure |3] Similarly, we 
can draw the face poset of Vg (see Figure |4] for 
the face poset of Vr 3 )- 

An n x n semi-magic square of magic sum 
r is an n x n matrix with nonnegative integer en- 



Figure 5: A magic labeling of T 3 and its corre- 
sponding symmetric magic square. 



Figure 6: Two different graph labelings associated 
to a semi-magic square. 



tries such that the entries of every row and column 
add to r. Doubly stochastic matrices are n x n 
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matrices in W 1 such that their rows and columns 
add to 1. The set of all n x n doubly stochastic 
matrices form a polytope B n , called the Birkhoff 
polytope. See [7 1, 1 8], or [ 15] for a detailed study 
of the Birkhoff polytope. 

A symmetric magic square is a semi-magic 
square that is also a symmetric matrix. Let H n {r) 
denote the number of symmetric magic squares 
of magic sum r (see 0, [19 1, and the references 
therein for the enumeration of symmetric magic 
squares). We define the polytope <S n ofnxn 
symmetric magic squares to be the convex hull 
of all real nonnegative n x n symmetric matrices 
such that the entries of each row (and therefore 
column) add to one. 

A one-to-one correspondence between sym- 
metric magic squares M = [rriij] of magic sum r, 
and magic labelings of the graph T n of the same 
magic sum r was established in fl9l : let ey de- 
note an edge between the vertex v.- t and the vertex 
Vj of T n . Label the edge of T n with rriij, then 
this labeling is a magic labeling of T n with magic 
sum r. See Figure |5] for an example. Therefore, 
we get Vr n = S n and H Tn (r) = H n {r). 

Corollary 1.4.1. The polytope of magic labelings 
of the complete general graph Vr n is an n(n — 
l)/2 dimensional rational polytope with the fol- 
lowing description 



L 

L(eij 



n(n + l) 
I 2 ; 



>0; 

1 < i,j < n,i < j, 

E}=1 L (ej;) + E"=i+i Hdj) 
fori = l,...,n 



The d-dimensional faces ofVr„ are d-dimensional 
polytopes of magic labelings of positive graphs 
with n vertices and at most n + d edges. There 
are C^™^ 1 ) faces of Vr 2n that are copies of the 
Birkhoff polytope B n . 



We define a digraph D to be a positive di- 
graph if the corresponding bipartite graph Gd is 
positive. 

Theorem 1.5. Let D be a positive digraph with at 
least one edge. Then, Vd is an integral polytope 
with dimension q — 2n + b, where q is the number 
of edges of D, n is the number of vertices, and b is 
the number of connected components of Go that 
are bipartite. The d-dimensional faces ofVr> are 
the d-dimensional polytopes of magic labelings of 
positive subdigraphs of D with n vertices and at 
most 2n — b + d edges. 

Let Il n denote the complete digraph with n 
vertices, i.e, there is an edge from each vertex to 
every other, including the vertex itself (thereby 
creating a loop at every vertex), then Gu n is the 
the complete bipartite graph K n ^ n . We get a one- 
to-one correspondence between semi-magic squares 
M = [rriij] of magic sum r and magic label- 
ings of Il n of the same magic sum r by labeling 
the edges of Il n with rriij. This also implies 
that there is a one-to-one correspondence between 
semi-magic squares and magic labelings of K n ^ n 
(this correspondence is also mentioned in [ 1 8 1 and 
l20l ). See Figure |6]for an example. 

A good description of the faces of Birkhoff 
polytope is not known [14|. We can now give an 
explicit description of the faces of the Birkhoff 
polytope. 
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Figure 7: The edge graph of the Birkhoff Polytope 
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Figure 9: The generators of the 2-dimensional 
faces of the Birkhoff Polytope B3. 



Figure 8: The generators of the edges of the 
Birkhoff Polytope B 3 . 



Corollary 1.5.1. Vu n is the Birkhoff polytope B n . 
The d-dimensional faces of B n are polytopes of 
magic labelings of positive digraphs with dimen- 
sion d, n vertices and at most 2n + d — 1 edges. 
The vertices ofPo, where D is a positive digraph, 
are permutation matrices. 

See Figure Qfor the edge graph of B3. Two 
faces of a polytope of magic labelings of a graph 
(or a digraph) are said to be isomorphic faces if 
the subgraphs (subdigraphs, respectively) defin- 
ing the faces are isomorphic. A set of faces is 
said to be a generating set of d-dimensional faces 
if every d-dimensional face is isomorphic to one 
of the faces in the set. See Figures [8] [TO] and 
[TTlfor the generators of the edges, the two dimen- 
sional faces, the facets, and the Birkhoff polytope 
B3, respectively. The numbers in square brackets 
in the figures indicate the number of faces in the 
isomorphism class of the given face. 

The proofs of Theorems 11.41 and 11.51 and 
Corollaries 11.4.11 and II .5. II are presented in Sec- 



[6] [3] 




Figure 10: The generators of the facets of the 
Birkhoff Polytope B 3 . 




Figure 11: The Birkhoff Polytope B3. 



tion|2] 

2 Polytopes of magic labelings. 

A polytope V is called rational if each vertex of 
V has rational coordinates. 

An element f3 in the semigroup Sc G is said 
to be completely fundamental, if for any positive 
integer n and a, a' G <5c*g> = a + a ' implies 
a = i/3 and a' = (n — i)j3, for some positive 
integer i, such that < i < n (see [ 17 1). 

Lemma 2.1. is a rational polytope. 

Proof. Proposition 4.6.10 of Chapter 4 in IT71 
states that the set of extreme rays of a cone and 
the set of completely fundamental solutions are 
identical. Proposition 2.7 in [18| states that ev- 
ery completely fundamental magic labeling of a 
graph G has magic sum 1 or 2. Thus, the extreme 
rays of the cone of magic labelings of a graph G 
are irreducible 2-matchings of G. We get a vertex 
of Vg by dividing the entries of a extreme ray by 
its magic sum. Thus, Vg is a rational polytope. □ 

Lemma 2.2. The dimension of Vg is q — n + b, 

where q is the number of edges ofG, n is the num- 
ber of vertices, and b is the number of connected 
components that are bipartite. 

Proof. Ehrhart's theorem states that the degree of 
Hg{t) is the dimension of Vg [6]. The degree of 
Hg{t) is q — n + b by Theorem ll.2l Therefore, 
the dimension of Vg is q — n + b. □ 

Lemma 2.3. The d-dimensional faces of Vg are 
the d-dimensional polytopes of magic labelings of 
positive subgraphs of G with n vertices and at 
most n — b + d edges. 

Proof. An edge e labeled with a zero in a magic 
labeling L of G does not contribute to the magic 
sum, therefore, we can consider L as a magic la- 
beling of a subgraph of G with the edge e deleted. 
Since a face of Vg is the set of magic labelings 
of G where some edges are always labeled zero, 
it follows that the face is also the set of all the 
magic labelings of a subgraph of G with these 



edges deleted. Similarly, every magic labeling 
of a subgraph H with n vertices corresponds to 
a magic labeling of G, where the missing edges 
of G in H are labeled with 0. Now, let H be a 
subgraph such that the edges e r \,. . . , e rm are la- 
beled zero for every magic labeling of H. Then 
the face defined by H is same as the face defined 
by the positive graph we get from H after delet- 
ing the edges e r %, . . . , e rm . Therefore, the faces 
of Vg are polytopes of magic labelings of posi- 
tive subgraphs. 

By Lemma l2~2l the dimension of Vg is q — 
n + b. Therefore, to get a d-dimensional polytope, 
we need to label at least q — n + b — d of G edges 
always 0. This implies that the d-dimensional face 
is the set of magic labelings of a positive subgraph 
of G with n vertices and at most n — b + d edges. 
□ 

The proof of Theorem l 1 .4l follows from Lem- 
mas 12.11 12.21 and 12.31 We can now prove Corol- 
lary OH 

Proof of Corollary \1.4.1\ 

It is clear from the one-to-one correspondence 
between magic labelings of T n and symmetric magic 
squares that Vr„ has the given description. Since 
the graph T n has n ( n 2 +1 ) edges and n vertices, and 
every graph is a subgraph of T n , it follows from 
Theorem ll.4l that the dimension of Vr n is n ( n ~ 1 ) ; 
the d-dimensional faces of Vr n are d-dimensional 
polytopes of magic labelings of positive graphs 
with n vertices and at most n + d edges. 

We can partition the vertices of T2 n into two 
equal sets A and B in ( 2n „ X ) ways: Fix the ver- 
tex v\ to be in the set A , then we can choose the 
n vertices for the set B in ( 2n ? 7 1 ) ways, and the 
remaining n — 1 vertices will belong to the set A. 
By adding the required edges, we get a complete 
bipartite graph for every such partition of the ver- 
tices of T2n- Thus, the number of subgraphs of 
T2n that are isomorphic to K n n is ( 2n ^ 1 ) ■ There- 
fore, there are { 2n ~ 1 ) faces of Vr 2 „ that are Birkhoff 
polytopes because every isomorphic copy of K nn 
contributes to a face of Vr 2n ■ d 

Lemma 2.4. Let G be a graph with n vertices. A 
labeling L of G with magic sum s can be lifted to 



a magic labeling L' of the complete general graph 
T n with magic sum s. 

Proof. Since G is a subgraph of T n , every labeling 
L of G can be lifted to a labeling V of Y n , where 



L'( eij ) 



L(eij) if eij is also an edge of G, 







otherwise. 



Since the edges with nonzero labels are the 
same for both L and L', it follows that the magic 
sums are also the same. □ 

Lemma 2.5. Let G be a graph with n vertices. 
The minimal Hilbert basis of Cq can be lifted to 
a subset of the minimal Hilbert basis ofCr n - 

Proof. If L is an irreducible magic labeling of G, 
then clearly it lifts to an irreducible magic labeling 
L' of T n . Since the minimal Hilbert basis is the 
set of all irreducible magic labelings, we get that 
the minimal Hilbert basis of Cq corresponds to a 
subset of the minimal Hilbert basis of Cr„ ■ □ 

Similarly, we can prove: 

Lemma 2.6. For a digraph D with n vertices, a 
magic labeling L with magic sum s can be lifted 
to a magic labeling V ofH n with the same magic 
sum s. The minimal Hilbert basis of Cd can be 
lifted to a subset of the minimal Hilbert basis of 
Cn„- 

Lemma 2.7. Let D be a digraph with n vertices. 
All the elements of the minimal Hilbert basis of 
Cd have magic sum 1. 

Proof. It is well-known that the minimal Hilbert 
basis of semi-magic squares are the permutation 
matrices (see [ 15 1) and therefore have magic sum 
1 . The one-to-one correspondence between magic 
labelings of II n and semi-magic squares implies 
that the minimal Hilbert basis elements of Cu n 
have magic sum 1. It follows by Lemma l2~6l that 
all the elements of the minimal Hilbert basis of 
C d have magic sum 1 . □ 

We now prove our results about polytope of 
magic digraphs. 

Proof of Theorem 17.51 By Lemma 12.71 all the 
elements of the Hilbert basis of Cd have magic 



sum 1. Since the extreme rays are a subset of the 
Hilbert basis elements, it follows that the vertices 
of Vd are integral. Since Vd = 'Pgb, it fol- 
lows by Theorem 11.41 that the dimension of Vd 
is q — 2n + b; the d-dimensional faces of Vd are 
the d-dimensional polytopes of magic labelings of 
positive subdigraphs of D with n vertices and at 
most 2n — b + d edges. □ 

We derive our results about the faces of the 
Birkhoff polytope as a consequence. 

Proof of Corollary \1.5.1\ The one-to-one 
correspondence between semi-magic squares and 
magic labelings of U n gives us that Vn n = B n . 
Since every digraph with n vertices is a subgraph 
of Il n , by Theorem ll.51 it follows that its d-dimensiona] 
faces are d-dimensional polytopes of magic label- 
ings of positive digraphs with n vertices and at 
most 2n — 1 + d edges. Since the vertex set of 
a face of B n is a subset of the vertex set of B n 
it follows that the vertices of Vd, where D is a 
positive digraph, are permutation matrices. □ 

3 Applications. 

In this section, we present some examples and ap- 
plications of magic graphs. Interesting examples 
of magic digraphs are Cayley digraphs of finite 
groups. Let G be a finite group {gi, g 2 , . . . , g n = 
I}. The Cayley group digraph of G is a graphical 
representation of G: every element gi of the group 
G corresponds to a vertex Vi (i = 1, 2, . . . , n) and 
every pair of distinct vertices Vi, Vj is joined by an 
edge labeled with a where g a = gjg^ 1 \ 12 1. For 
example, the Cayley digraph for the permutation 
group 

5 3 = {5i = (123), 52 = (132), 93 = (23), 

54 = (12), 55 = (13), 56 = /} 

is given in Figure [T2l 

Proposition 3.1. The Cayley digraph of a group 
of order n is a magic digraph with magic sum 

n(n—l) 
2 • 

Proof. Let denote an edge between the vertex 
Vi and Vj of the Cayley digraph such that V{ is the 



Figure 13: Graph Factorization. 



Figure 12: Cayley digraph of the group S3 [ 12 1. 



initial vertex and vj is the terminal vertex. Let v\ 
be a vertex of the Cayley digraph, and let a be an 
integer in the set {1, 2, . . . , n — 1}. Let g p = g a gi 
and let g q = gig a . Then, the edges e/ p and e q \ are 
labeled by a. Also, gjg^ 1 = g n = I if and only if 
i = j. Hence, a Cayley group digraph is a magic 

digraph with magic sum 1 + 2 + \- (n — 1) = 

V^Hzll ( see also Chapter 8, Section 5 in 02 D- □ 

A digraph is called Eulerian if for each ver- 
tex v the indegree and the outdegree of v is the 
same. Therefore, Eulerian digraphs can also be 
studied as magic digraphs where all the edges are 
labeled by 1 (see [4| for the applications of Eule- 
rian digraphs to digraph colorings). An n-matching 
of G is a magic labeling of G with magic sum at 
most n and the labels are from the set {0, 1, ... , n} 
(see lTi3l . chapter 6). A perfect matching of G is 
a 1 -matching of G with magic sum 1. 

Proposition 3.2. The perfect matchings of G are 
the minimal Hilbert basis elements of Co of magic 
sum 1 and the number of perfect matchings of G 
isH G (l). 

Proof. Magic labelings of magic sum 1 always be- 
long to the minimal Hilbert basis because they are 
irreducible. Therefore, perfect matchings belong 
to the minimal Hilbert basis because they have 
magic sum 1. Conversely, every magic labeling 
of magic sum 1 is a perfect matching. So we con- 
clude that the perfect matchings of G are the min- 



imal Hilbert basis elements of Cq of magic sum 
1 . The fact that the number of perfect matchings 
of G is Hg{1) follows by the definition of Hq(1). 
□ 

Hilbert basis can also be used to study fac- 
torizations of labeled graphs. We define Factors 
of a graph G with a labeling L to be labelings 
Li,i = 1,..., r of G such that L(G) = J2i=i L i( G )> 
and if Lj(efc) 7^ for some edge of G, then 
Lj(ek) = for all j ^ i. A decomposition of L 
into factors is called a. factorization of G. An ex- 
ample of a graph factorization is given in Figure 
□J See Chapters 1 1 and 12 of [ 12] for a detailed 
study of graph factorizations. 

Our results enable us to reprove some known 
facts about the Birkhoff poly tope as well. For 
example, Theorem 11.51 gives us that the dimen- 
sion of B n is (n — l) 2 . The leading coefficient 
of the Ehrhart polynomial of B n is the volume 
of B n . This number has been computed for n = 
1, 2,..., 9 (see 15] and d). 

The software 4ti2 ATI can be used to find 
the Hilbert bases of Cq, and the software LattE 
fTOll can be used to compute the generating func- 
tions of Hc{r) effectively. See [3| for other re- 
sults about magic labelings of graphs. See 1221 
for a study of other types of magic graphs. 
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